Light scattering by optically anisotropic scatterers I: 
T— matrix theory for radial and uniform anisotropies 

A.D. Kiselev/'Q V.Yu. Reshetnyak,2'| and T.J. Sluckin^'l 

^Chernigov State Technological University, 
Shevchenko Street 95, 14027 Chernigov, Ukraine 
^Kiev University, Prospect Glushkova 6, 03680 Kiev, Ukraine 
^Faculty of Mathematical Studies, University of Southampton, Southampton, SOU IB J, UK 

(Dated: February 2, 2008) 

Abstract 

We extend the T-matrix approach to light scattering by spherical particles to some simple 
cases in which the scatterers are optically anisotropic. Specifically we consider cases in which the 
spherical particles include radially and uniformly anisotropic layers. We find that in both cases the 
T-matrix theory can be formulated using a modified T-matrix ansatz with suitably defined modes. 
In a uniformly anisotropic medium we derive these modes by relating the wave packet representation 
and expansions of electromagnetic field over spherical harmonics. The resulting wave functions are 
deformed spherical harmonics that represent solutions of the Maxwell equations. We use these 
modes to express the equations for the T-matrix elements in terms of computationally tractable 
coefficient functions. 
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I. INTRODUCTION 



The problem of light scattering by particles of one medium embedded in another has 
a long history, dating back almost a century to the classical exact solution due to Mie [|l|. 
The Mie solution applies to scattering by uniform spherical particles with isotropic dielectric 
properties. More recently this strategy has been successfully applied to ellipsoidal particles 
and some circumstances in which the dielectric tensor is anisotropic 0, ^, ||, ||, 0, |^. 

Unfortunately, Mie-type solutions, although exact, are not always physically meaningful. 
An alternative strategy useful in the long wavelength limit is the so-called Rayleigh-Gans (R- 
G) approximation 0], which is closely allied to the classic Born approximation of quantum 
mechanics. However, detailed information on the precise range of validity of this approxi- 
mation often requires solution of the Mie problem, which is only available in some specific 
cases. 

There are a large number of physical contexts in which it is useful to understand light 



scattering by impurities ||T0[- A particular example of recent interest concerns liquid crystal 
devices. There are now a number of systems in which liquid crystal droplets are suspended 
in a polymer matrix - the so-called PDLC systems - or the inverse system, involving colloids 
now with a nematic liquid crystal solvent. These inverse systems are commonly known as 
filled nematics 0, 0. 



In such systems one needs to calculate light scattering by composite anisotropic particles 
embedded in an isotropic or an anisotropic matrix. The model scatterer usually consists of 
a small central isotropic particle ("the core"), coated by a much larger region in which the 
optical tensor is anisotropic. This is equivalent to examining light scattering by a composite 
particle consisting of the central core plus a surrounding liquid crystalline layer. 

The analysis of a Mie-type theory uses a systematic expansion of the electromagnetic 
field over vector spherical harmonics. The specific form of the expansions is known as the 
T-matrix ansatz. This has been widely used in the related problem of light scattering by 
nonspherical particles [113], n . 



Recently in [|^ |^ we have studied the scattering problem for the optical axis distributions 
of the form: n^f -|- n^'d + n^p^p. By using separation of variables and expansions over vector 
spherical harmonics, we have developed the generalised Mie theory as an extension of the T- 
matrix ansatz [1^. This theory combines computational efficiency of the T-matrix approach 
and well defined transformation properties of the spherical harmonics under rotations. 

In this paper we discuss this theory in more detail and explain how this approach can be 
extended to the case of uniformly anisotropic spherical particles. 

The layout of the paper is as follows. General discussion of the model is given in Sec. II. 
Then in Sec. fTT] we outline the T-matrix formalism for the isotropic medium in the form 
suitable for subsequent generalisation. In Sec. |IV|, as the simplest case to start from, we 
consider how the T-matrix ansatz applies for the radially anisotropic layer. We find that 
the structure of electromagnetic modes in the layer requires modification of the standard 
T-matrix ansatz. In addition, we detail computing the elements of T-matrix. 

In Sec. we describe the method to put the scattering problem into the language of 
T-matrix by linking the representations of plane wave packets and of spherical harmonics. 
For uniformly anisotropic scatterer we define generalised spherical harmonics and show that 
the effect of angular momentum mixing can be treated efficiently. 

Finally in Sec. ^ we draw together the results and make some concluding remarks. In 
particular, we emphasise the importance of anisotropy effects by making comparison between 
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angular distributions of scattered wave intensities for radially anisotropic layer and effective 
isotropic layer of the same scattering efficiency. 

Details on some technical results are relegated to Appendices ^ - |C[ 

II. MODEL 

We consider scattering by a spherical particle of radius Ri embedded in a uniform isotropic 
dielectric medium with dielectric constant eij = eSij and magnetic permeability fiij = fiSij. 
The scattering particle consists of an inner isotropic core of radius R2, surrounded by an 
anisotropic annular layer of thickness d = Ri — R2. 

Within the inner core of the scatterer the dielectric tensor e, and the magnetic permittiv- 
ity /X take the values Sij = e2Sij, fi^j = iJ,2^ij. The dielectric tensor within the annular layer 
is locally uniaxial. The optical axis distribution is defined by the vector field n. (Hats will 
denote unit vectors.) Then within the annular layer ejj(r) = ei5ij + Aei(n(r) ® n(r))jj and 
Mij = A^i^ij- The unit vector n corresponds to a liquid crystal director for material within 
the annular region. 




Isotropic core 



FIG. 1: Distributions of optical axis in the anisotropic layer around a spherical particle for radial 
and uniform structures: a) n = f ; b) n = z 

We shall suppose that the director field can be written in one of the following forms 

n = z (la) 
fi = f (lb) 

where f = (sin 9 cos 0, sin 9 sin 0, cos 9) is the unit radial vector; and 9 are Euler angles of 
the unit vector f . x, y and z are the unit vectors directed along the corresponding coordinate 
axes. 

In Fig. m we have shown these director distributions. Fig. Qa shows the radial (and 
spherically symmetric) director distribution. The scattering problem in this case has already 
been discussed by Roth and Digman Q], and in this case the spherical symmetry of the 
problem plays an important role in rendering the Maxwell equations soluble. 

We show in Fig. [l|b the case in which the optical axis is directed along the 2;-axis and 
is uniformly distributed within the annular layer; this is given by Eq. (p^. The case where 
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the scatterer is a long cylinder parallel to fi presents no difficulties and can be treated in 
cylindrical coordinate system [1^. Scattering from spherical uniformly anisotropic particles 



is not exactly soluble and has been studied by using the Rayleigh-Gans method and 



the anomalous diffraction approximation in |T^, |T1 

A simple limit of the physical situations we consider puts e = ei = €2, with ji = ni = 
fi2- The first condition allows us to concentrate on situations in which the scattering is 
governed by the anisotropic part of the dielectric properties. This distinguishes our case 
from other studies of scattering by spheres, in which the isotropic optical contrast dominates. 
However, there is also a motivation for this hypothesis in terms of liquid crystal device 
physics, and we shall discuss this at greater length in a subsequent paper. However, the 
result of the hypothesis is that the scattering by our model spheres disappears in the limit 
of zero anisotropy. 



III. r MATRIX APPROACH IN ISOTROPIC MEDIUM 
A. T matrix ansatz 

In this subsection we remind the reader about the relationship between Maxwell's equa- 
tions in the region of a scatterer and the formulation of scattering properties in terms of the 
T-matrix |TU[. Our formulation is slightly non-standard. Some technical details, which 
can be omitted at ffist reading, have been relegated to the appendices. 

We shall need to write the Maxwell equations for a harmonic electromagnetic wave (time- 
dependent factor is exp{—iujt}) in the form: 

-iriilfiiki]'^ V X E = H , (2a) 
ifiilriiki]'^ V X H = E + Ui{h ■ E)h . (2b) 

where rij = y/eijli are refractive indexes for the regions, where i?2 < ''^ < -Ri (^ = 1) and 
r < R2 {i = 2); ki = Uik^ac {k^ac = uj/c = 27i/X is the free-space wavenumber). We define 
the anisotropy parameter as Ui = Aei/ei (in the annular layer). Then inside the isotropic 
core U2 = 0. Finally, for brevity, in the region outside the scatterer r > Ri, the index will 
be suppressed, giving k = nk^ac and u = 0. 

The electromagnetic field can always be expanded using the vector spherical harmonic 
basis, Yj+sj7n{4>,6) = Yj+sjm{r) [6 = 0,±1) |T8[, as follows: 

E = J2e,^ = J2 [pfi(OYj°i(f) +p(l(r)Y(;l(f) +p5:V)Yj^ , (3a) 

jm jm 

H = 5^ H,„ = 5^ [C(r)Yfi(f) + q^^)Y^ir) + qlZ\r)Y\Z\f)] , (3b) 



are 



where y\Z^ = Y„-^ and Y^ = [j/{2j + 1)]1/2Y,+i,„ + [{j + l)/(2j + Ip'Y,^,,^ 

electric and magnetic harmonics respectively, and Yj^ = [j/{2j + l)]"'^/^Yj_i — [{j + 
l)/(2j -|- l)]^/^Yj_|_i are longitudinal harmonics (a number of relations for the vector 
spherical harmonics used throughout this paper are considered in Appendix The electric 
field is now completely described by the coefficients {p^j^{r)} and similarly the magnetic field 
is now described by {qj^if')} with a = {o, e, m}. 
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In order to find the coefficient functions we use separation of variables. This imphes 
that the expansions (j^) must be inserted into Maxwell's equations (0). The coefficient 
functions then can be derived by solving the resulting system of equations. In the simplest 
case of an isotropic medium the coefficient functions can be expressed in terms of spherical 



Bessel functions, jj{x) = [n / {2x)Y^'^ Jj^i/2{x), and spherical Hankel functions [T^, hj 



XI 



['k/{2x)Y^'^Hj_^-^^2{x), and their derivatives as follows 

E,^ = a.^Mj.;::) (p, f ) + ^,^M^^ (p, f ) - ^ («,n.M(.S,(p, f ) + 4-„m51(p, f )) , (4a) 
H,^ = a,-„M('::)(p, f) + P,mM^;:^ip, f) + ^ [a^mM^^XP, f) + P,mM^;iip, f)) , (4b) 

where ajm, ajm, (3jm and are integration constants; the vector functions M^-^ and M^-^ 
are given by 

Mt^(p,f)=j,(p)Yj:)(f), 

M5.:i(p, f ) = /^j,(p) Yfi(f) + V^2IZ±iI J- .(p)Yj°2(f ) , (5a) 

M^iP. r) = Dhf\p)Y\'iir) + V^II+lI hf\p)Y\t{r) , (5b) 

where Df(x) = x~^ — (xf{x)) and p = kr. 
ax 



There are two cases of Eq. (^aj) that are of particular interest. They correspond to the 
incoming incident wave, {EincHinc}, and the outgoing scattered wave, {EscayHsca}'- 

jm jm jm vr' / ^ jm jmxr') J ' 

H(mc) ~ (mc) -« «-(m) / ^\ . ^ (inc) /f(e) / ^\ 
jm = (^jm M} J (p, r) + - a}„ M),^{p, r) , (6) 

P 

^ = (p, f ) - ^ ^]r^ m(1(p, f ) . 

H$r^ = /3]r^ iP, + l & m(1(p, f ) . (7) 

p 

Now the incoming incident wave is characterised by amplitudes a^j^"* ■> '^fm^'' and the scattered 
outgoing waves are similarly characterised by amplitudes p'f^\ Pjm"^^- Our task is now to 
relate {a} and {/?}. 

In this regime, a transverse plane wave incident in the direction is specified by an unit 
vector kinc, with 

E,„e = E^^"'^) exp(^ kinc ■ r) , E^^"^) = 4'"'^e,(k,„e) , k,„e = kk^^c • (8) 
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We show in ( |B5| ) the coefficients {a} of the expansion @ takes the form: 



u=±l 



~ (inc) 



z/=±l 



where -D^^, is the Wigner D-function and the basis vectors e±i(kj„c) are perpen- 

dicular to kj„c and defined by Eq. (p2D. 



Thus outside the scatterer the electromagnetic field is a sum of the transverse plane wave 

){inc) n(inc 

jm r'jm 



incident in the direction specified by an unit vector kj„c iPjm'^^ = Pfm''^ = 0) and the outgoing 



wave with aj^""* = af^^ = as required by the Sommerfeld radiation condition. 

So long as the scattering problem is hnear, the coefficients /S]^"'* and /S]^"'* can be written 
as linear combinations of a]*^'^'' and a^j^^'- 



p. 



(sea) 
jm 



\_ jm,j'm' ^j'm' ' ^ ''-jm,j'm' ^j'm' 



j,m,' 



J' ,m' 



'^rp2l Hnc) rp22 ~{inc) 

jm, j'm' ^j'm' ' jm, j'm' ^j'm' 



(10) 



These formulae define the elements of the T-matrix in the most general case. 

In general, the outgoing wave with angular momentum index j arises from ingoing waves 
of all other indices j'. In such cases we say that the scattering process mixes angular 
momenta fl^. The light scattering from the uniformly anisotropic scatterer, depicted in 
Fig. |l]b, provides an example of such a scattering process. In this case the cylindrical 
symmetry of the optical axis distribution causes the T-matrix to be diagonal over azimuthal 
indices m and m': 
using matrix notations: 



nn 

jiT^t j'm' 



iJmm' -'-jj'-m- 



Then we can conveniently rewrite the relation (|ToD 



I ^Ic^ca) = E ^'^jj'-.rn , ^'%nc) > [^\jj';m = T^'l"" ^^2'"' • (H) 
f^/nl3}.^'J \f^/n(^j'mj Vn';m -^jj';mj 

In simpler scattering processes, by contrast, such angular momentum mixing does not 
take place. Many quantum scattering processes and classical Mie scattering belong to this 
category. It is seen from Fig. Qa that radial anisotropy keeps intact spherical symmetry of 
the scatterer. The radially anisotropic annular layer thus exemplifies a scatterer that does 
not mix angular momenta. The T-matrix of a spherically symmetric scatterer is diagonal 
over the angular momenta and the azimuthal numbers: Tj^"^^' — ^jj'^mm'Tj^^ . 



B. Scattering Amplitude Matrix 

We have seen the relation between the scattered wave (0) and the incident plane wave (P) 
is linear. In the far field region (p ^ 1), where the asymptotic behaviour of the spherical 
Bessel and Hankel functions is known [p!9|] : 
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V+^hf{p),VDhf{p) ~ exp(«p)/p. 



(12) 



The scattering amplitude matrix A(ksca, kj„c), which relates E^ca and the polarisation vector 
of the incident wave E*^*"'^-' is defined in the following way 0, [TO, [TH]: 



4^'^'^) = (e;(k,,J, E,,J = p-^ exp(«p) 5^ A,,,(k,,„ k,^,)^!"'^) , z/ = ±1 (13) 



where an asterisk indicates complex conjugation, k^ca = f and e±i(ksca) = =F(''^ ±«<^)/"\/2- 
Eqs. (^, (|12D and the vector spherical harmonic relations Eqs. ([A9| ) can now be combined 
to yield the expression for the scattering amplitude matrix in terms of the T-matrix: 



KAKca. knc) = E E + l)(2j' + l)]^/'/^^;(k.ea)I?£v'(k.nc) " 

jm j'm' 

■luu'T^^ -iuT^'^ -\-7u'T'^^ ,7^22 1 (-iA\ 

[ jm,j'm' ^ jm,j'Tn' ' ^ jm,j'm' ' \ ' \ J 

For a spherically symmetric scatterer, this result can be expressed in the simplified form: 



^uu' {}^sca^ ^inc) ^ ^ -^jyjy' (k^ca; ^iric) 
3 

3 

(15a) 

DiA^sca, knc) = E ^^:(ksca)I?L,(k»c) • (15b) 

Eq. ( [L5b ) shows that the scattering amplitude matrix (|15a|) depends only on the angle 
between kj„c and kj^ic- 

All scattering properties of the system can be computed from the elements of the scatter- 
ing amplitude matrix. In Eq. (|1^ we see that this can be defined in terms of the elements of 
the T-matrix defined in Eq. ([T0|). Thus computation of these matrix elements is of crucial 
importance. 



C. Scattering Efficiency 

In order to find the total scattering cross section Csca we need to calculate the fiux of 
Poynting vector of the scattered wave S^'*'^"^ = c/ (Svr) Re(Esca x H *^^) through a sphere of 
sufficiently large radius and divide the result by IS'-*"'^-'!. 

From Eqs. (^ and (0) the asymptotic behaviour of the scattered outgoing wave in the 
far-field region is given by 
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2j + l 



in 



1/2 



/^jr^ - ^/^^ '/^jr^ ) ^m*-i(ksca) e_i(k 



-1 o{sca) 



(16a) 



.j+ljj(.ca) _ 



2j + l 



57r 



/5]r^ ) ^^:-i(k.ea)e_i(k 



/^jr^ - ^ ^^/^jr^ ) ^m* 1 (ksca) e+i (k. 



-1 /o{sca) 



(16b) 



where the relations between the vector spherical functions and Wigner D-functions (|A9|) 
have been used. 

The asymptotic relation ( |16a| ) give the relations (p!3|) and (|1^) that define the scattering 
matrix A. Similarly, from Eq. ( |16b|) we have the following relation for magnetic field of the 
scattered wave in the far field region: 



(mc) 



u'=±l 



(17) 



From Eqs. (p!3|), (p!7|) and ( |A3| ) we can now readily express the Poynting vector in terms 
of the scattering matrix: 



ly, u'=±l 



(18) 



where the superscript f indicates Hermitian conjugation. 

Using these expressions and the orthogonality relation ( Aid ), we can immediately deduce 
the result for the total scattering cross-section: 



^sca — rv / ^ \yjm I ^ \H'"' Hjm \ 



(19) 



where line = \Ej^ 



(mc) I 2 



u=±l 



Integrating the product of matrices that enter Eq. ( pISD over the angles of scattered wave 
gives 



[A ■ At] = 2nJ2Yl + 1)(2J' + 1)]'/' ^^Uk^c) D^Jk^nc) ■ 

jm j'm' 

[{luT'' + Ti2)t . {^v'T'' + T^2) + (^z.T^^ + T^')^ ■ {lu' T^' + T^')]^^^^,^, , (20) 
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where {f)sca= / d(f)sca / sin 6 sea d6 sea f- 

Jo Jo 

So, we have the matrix (A ■ A'^)sca which is, in general, non-diagonal and depend on 
incident wave angles. For a spherically symmetric scatterer this matrix is diagonal and 
independent of kj^c- In this case the cross-section (|I9|) can be written in the following form: 



Csea = ^Tr{A-A^)sca. (21) 

Note that for unpolarised incident light this relation holds even if a scatterer is not spherically 
symmetric. 

We can now relate the scattering cross-section and the elements of the T-matrix. For 
brevity, we restrict ourselves to the spherically symmetric case, considered at the end of 
the previous subsection. More precisely, we consider the scattering efficiency, Q, that is the 
ratio of Cgca and area of the composite particle, S = vri?^. For diagonal T-matrix, Eqs. ( pTD 
and ( ppp give the known result [0: 

^ 9 2 

Q = ^ = ^E E(2i + i)l^rr- (22) 

1 j=l m,n=l 

In order to characterise angular distribution of scattered light intensity let us suppose 
that the incident wave is linearly polarised and is propagating along the z-axis, kj„c = z. 
The wave vectors k^c and k^ca define the scattering plane. From Eqs. (|l^) and (|l^) we can 
find the amplitudes of the scattered wave components that are parallel {Ejf^"''') and normal 
(^j^{sea)^ to the scattering plane: 

W'^r^^Osea) = hnel\\{9sea) COs' ^ , (23a) 
(^-") = I^ne^±{Osea) sin^ ^ , (23b) 

i\\ = \ii + , 2_L = |ii — , (23c) 

Z,{9sea) = {kRi)'' + 1/2) dUOsea) [T}' + uTf], (23d) 

i 

where ip is the angle between the polarisation vector of the incident wave and the scattering 
plane. 

In subsequent sections we shall use the intensity iseai^sea) 

iseaiOsea) = {^\\{0sea) + i±{Osea)) /2 (24) 

and the factor characterising the degree of depolarisation (depolarisation factor) PdepiGsea) 

Pdep{Osea) = 1 - P{6sea), P{6sea) = , (25) 

iWlosea) + i±['ysea) 



as quantities characterising angular distribution and polarisation of the scattered wave I^T 
Note that averaging iseaiOsea) over the scattering angle gives the scattering efficiency: 

/*7r 

Q / '^seai^sca) SIB. Ogea d.9sca ■ 
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IV. SCATTERING FROM RADIALLY ANISOTROPIC LAYER 



In Sec. [Ill A| we started from the general expansion for electromagnetic field over the 
vector spherical harmonics (^). Then the fields in isotropic medium were expressed in terms 
of the modes, M^-^-* and M^-^^ (see Eq. @). This expression is known as the T-matrix 



ansatz O, 14 



We shall write down the results for electromagnetic field within the radially anisotropic 
layer as they are given in [Q. These can be written in the form similar to the T-matrix 
ansatz: 



Ej^ = a,^vf2 (P, i) + /3,mP5™^ (P, f ) - ^ («,,nPfi(p, f ) + km^fAP. f )) > (26a) 

= «,,.Q;:^ (p, f ) + /^..^qJ™) (p, f ) + ^ («,,nQ51(p, f ) + Mfiip, f )) . (26b) 

For radial anisotropy the modes that enter Eqs. (|26| ) are given by 

= rM) ' pS = ^Jj(Pi) YS(f ) + [jG + Pi' jM) Y]l(f ) , (27a) 

Pt^ = M^\p^,^), Q5-l = M(.l(pi,f), j(j + l)=j(j + l)/(l + ^i)- (27b) 

In the next section we shall find that the fields in uniformly anisotropic medium can also 
be written in the form (^). However, the expressions for the normal modes in this case 
will differ from those of Eqs. (p7|). In both cases the modes represent solutions of Maxwell's 
equations and turn into the isotropic medium modes (^) in the limit of small anisotropy 
parameter, u ^ 0. 

The radial anisotropy harmonics ( ]26| ) do not involve angular momentum mixing and as 
such they only have contributions proportional to vector spherical harmonics with angular 
momentum numbers given j and m. The important difference between the two cases lies in 
the fact that for the uniformly anisotropic layer this is no longer true. 

The fields inside the isotropic core of the particle similarly also involve no angular mo- 
mentum mixing: 



^ = «S mJ-) (p„ f ) - g a]l M^{p„ f ) , (28a) 
= «jl M^{p„ f) + ^ cy^ M;.1(p2, f) , (28b) 

where p, = kiV. 

We now pass on to the calculation of the scattering cross-section using the T-matrix 
method. The formulae presented in this subsection form a key element in the input to this 
calculation. 



A. T matrix: radial anisotropy 

In order to calculate the elements of T-matrix, we need to use continuity of the tangential 
components of the electric and magnetic fields as boundary conditions at r = i?2 and r = Ri. 
Equivalents, the functions pj.^] = (Yj^^* ■ E)^, and gfi = (y]^^* ■ H)f, a = m, e, must be 
continuous at each boundary: 
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i;JliRi + 0) = p]''J,iR, - 0), + 0) = qf^{R, - 0), z = 1, 2 ; a = m, e . (29) 

For the radially anisotropic distribution (|l^) different angular momentum are decoupled. 
We have shown elsewhere that the T-matrix can be computed in the closed form and 
here we give more details of this calculation. 

Expressions for {p^^, q^^} ambient medium have been given implicitly in Eqs. 

and (|^). Similar expressions in the annular layer have been given by Eqs. (|2^) and (P7\). 
Finally, the analogous expressions inside the isotropic core have been given in Eq. (p^). 

In order to determine the T-matrix we shall insert these expressions into the boundary 
conditions (^Qf ). This will yield a system of eight linear algebraic equations for the ten 
quantities: af^''\ a^-Z^\ (3'f^"'\ a^^, (3jm, "im, Pjm, afi and af^. After eliminating 

all internal variables we shall be left with two equations in the four unknowns: aj'm \ oi"jm\ 

l^frn^i i^fm^ ■ These equations will define the T-matrix. 

In order to carry out this procedure efficiently we combine Eqs. and ( p6D by using 
the compact matrix notation for the components inside the anisotropic layer: 



where 



Rfr) 



(m) I 



R(r) 



(3jn 



a 



jm 



ni/ii ^Djj{pi) ni/ii ^Dh) '{pi) 




V 














jjipl) 

-fiin^^Dj-Api] 








h^\p^ 



\ 



-pin^^Dh'j\pi)j 



(30) 



(31) 



We now apply the boundary conditions at the dielectric discontinuities on the inside and 
outside of the anisotropic layer. In the matrix notation this yields: 



R2 



Ri 



'jm 



a 



( bi(P2)]2 \ 



a 



(c) 
jm 



n2/i2Mjj(P2)]2 






0^- 



(sea) 
jm 



I 

I \hf{p)V \ 

nlp\hf{^py{, 


V / 



+ 3^- 



, (mc) 

+ a- 



I Up)]i \ 

n/p[jj{p)]'i 





. ~ (mc) 



/ \ 



_bi(P2)]2 
V2[ji(p2)]2y 

/ ° 



(sea) 

[hf{p)] 
\-p/n hf{p)^J 

( \ 



bi(p)]i 



+ 



(32a) 



(32b) 
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where Rj = R\r=R,, Df{x)\r=R, = [/(a;)] - and f{x 
The amphtude in the anisotropic layer, {ajm, P 



,Pjm}, can be related to the am- 



plitudes inside the core, {ctj^m, ctfri} , by multiplying both sides of Eq. ( |32aD by R2^. Sub- 



stituting this result into the left hand side of Eq. ( |32b|) enables the amplitudes inside the 
anisotropic layer to be eliminated, resulting now in a system of four equations for the six 



amplitudes {a 



(c) ~ (c) 

■ (y ■ 



n{sca) M 
r'jm 5 r'j 



(sea) 
m 



a 



(inc) ~ (inc) 



jm 



a 



jm 



where 





C12 


-[hf{p)]i 


\ 






C21 


C22 ~ 









~(c) 
01- 
jm 


C31 


C32 





-[hf{p)]i 




n{sca) 


\C41 


C42 





f^/n[hf\p)]'J 








/ 


[Ji(p)]i \ 


( 








(inc) 

a ■ 

jm 



^//^bj(p)]i 






. ~ (inc) 

-\- a ■ 





[Ji(p)]l 

\-f^/n[jjip)]'lJ 



c 



/Cll Ci2\ 

C21 C22 

C31 C32 

\C41 C42/ 

7m? ^jmJ 



RiR,2 



/ bi(P2)]2 
n2P2^[jj{P2)]2 



\ 



\ 



_[ji(P2)]2 
-"-2 V2[ji(P2)]2y 



(33) 



(34) 



We now eliminate {aj^, a^^} in the system (^). This yields the standard general form 
for the relation between the amplitudes of the incoming and outgoing waves (p^Of ), with the 
following explicit expressions for the elements of the T-matrix: 



yii ^ C2[jjip)]i-n/fici[jj{p)][ ^22 ^ C2 [jj{p)]i + p/ndi [jj{p)]i .g^s 
' n/pc,[h'f\p)][-C2[h'~;\p)],' ' p/nc^[hf\p)][ + ~C2[hf\p)l' 



where d = cn and q = Ci+2,2- 

We note that in this result the off-diagonal elements of the T-matrix vanish: T, 
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Tj^ = 0. This result is quite general for spherically symmetric scatterers, and physically 



means that there is no coupling between transverse electric and transverse magnetic waves. 
Equivalently there is no depolarisation scattering in this scatterer. 



Mathematically the result follows because the matrix (31) is block diagonal. Thus from 



Eq. (1^) we have C31 = C41 = C12 = C22 = 0, from which the result follows. The T- 
matrix is diagonal over the angular momenta j and the azimuthal numbers m: T^y^rnm' — 

S S A Tnn 
^jj'^mm'^nn' j 



V. SCATTERING FROM UNIFORMLY ANISOTROPIC LAYER 

By contrast with the case of the radially anisotropic scatterer considered in the previous 
section, the electromagnetic field inside the uniformly anisotropic annular layer can easily 
be described in terms of plane waves. Unfortunately this does not render the scattering 
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problem soluble. The difficulty lies in satisfying the boundary conditions for a spherical 
scatterer using the plane wave solutions of the Maxwell equations. 

The starting point of the T-matrix approach to this case involves finding a representation 
of the electromagnetic wave analogous to the generalised T-matrix ansatz (pGl). However, 
for symmetry reasons, the structure of the modes representing the fields in the uniformly 
anisotropic medium differs from that in the spherically symmetric isotropic (or a radially 
anisotropic) geometry. In particular, the lack of spherical symmetry implies that the angu- 
lar momentum number j is no longer a good quantum number. However, the cylindrical 
symmetry still guarantees conservation of the azimuthal number m. 

The procedure is as follows. In Sec. |V A| we provide methods of defining modes in a 
uniformly anisotropic material. These modes are: (a) solutions of the Maxwell's equations 
and (b) deformations of the isotropic spherical harmonics. The latter condition means that 
the isotropic modes Eqs. (§) are recovered in the weak anisotropy limit, u —>■ 0. Then 
in Sec. |VB| , we shall use these "quasi-spherical" wave functions to derive equations which 



enable the elements of T-matrix to be computed for a uniformly anisotropic annular layer. 



A. Angular Momentum Representation in the Anisotropic Layer 

We have seen above that solutions to Maxwell's equations can be written in two ways. 
Either they can be expressed as plane waves, or, using a separation of variables approach, 
they can be written as expansions over spherical harmonics. Deriving an expression for the 
T-matrix will require us to make a connection between these two alternative expansions. 
In this subsection we carry out this task. 

To do this we begin with an isotropic medium. In this special case both the spherical 
harmonics and the plane wave solutions are known. The result is a relation between the 
plane wave packets and the spherical harmonics. The procedure will then be generalised to 
cover the case of a uniformly anisotropic medium, so as to derive a set of " quasi-spherical" 
normal modes. 



1. Spherical modes and plane waves in isotropic media 

We start with the Maxwell equations (0) for an isotropic medium. Wave-like solutions to 
these equations written in terms of spherical coordinate basis functions are given by Eqs. 
and (P). An electromagnetic wave can alternatively be written as a superposition of plane 
waves: 



E = (exp(2pk ■ f) [E,{k) e,(k) + Ey{k) e^(k)])fc , (36a) 



H = - (exp(^pk ■ f) [E,(k) ey{k) - Ey{k) e,(k)] , (36b) 



/•27r PTT 

where (/)£= / d0fc / sin6'fcd6'fe/. 
^0 Jo 



'0 Jo 

In Appendix B we derive expressions connecting the isotropic spherical modes (|5a]) and 



the vector plane waves occurring in the superposition (|36D 
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= i"^'(47r)"V7r(2j + 1) (exp(ipk ■ f) 
M^ip, f ) = 



'k ' 



where 



D^-^(k) = exp(-2m0fc) (i' 



jm\ 



D]^i(k) ^ exp(-^m0,) d^J,{e,) = Dl^_,{k) + Dl^.ik) , D^(k) ^ Z}^,o(k) • 



Ay) 



(37a) 
(37b) 



(38a) 
(38b) 



These relations can be exphcitly verified by substituting the expansions ( [B6D into the right 
hand sides of Eqs. (p^). The hnear combinations of the modes Mj™''(p, f) and Mj^(p, f ) 
which enter the electromagnetic field harmonics (^) can now be expressed as a superposition 
of plane waves: 



a.^Mj^ (p, f) - (^a,M;i{p, f) = (exp(zpk ■ f) [E^'2i^) e,.(k) + E^^ik) e,(k)] , (39a) 
«.n.M(::V, f ) + -«.™M(.:i(p, f ) = - (exp(.p k . f ) [EfJ(k) e,(k) - E^i{k) e.(k)] , 

(39b) 



where 

i^]::J(k) = ^ [vr(2j + 1)]V2 {a,^D(^2*(k) - ^ ^ a.^^fj^lk)} , (40a) 
E^tik) = [vr(2j + 1)]V2 {za,^^Dfj*(k) + ^ a,„.I)fi*(k)} . (40b) 

We now sum Eqs. ( ^9| ) over j and m. This enables the amplitudes Ex{k) and Ey{k) in 
Eqs. (36) to be expressed in terms of Wigner D-functions: 



E^{k) = J2E\i.i^), Ey{k) = J2E\i{k). (41) 

jm jm 

This procedure has started from plane waves ( PB| ) and spherical harmonics (^), and finished 
with Eqs. (|4(]|)- (|4T|). This equation defines a basis set in the space of the angular dependent 
amplitudes. 

In fact we shall need to carry out the inverse process. The inverse process uses the 
expansions (|41|) to derive the expressions for the spherical modes M^-^ and Mj-^ from 
superpositions of plane waves ([361) . The procedure works as follows: 

(a) We substitute the expansions of the amplitudes -^^^(k) and Ey{k) from Eqs. (|4T|) into 
the superpositions (^). 
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(b) From the expressions for the electric (magnetic) fields we obtain the spherical modes 
as coefficient functions proportional to ajm and —fi/najm {(ijm and n/ fiajm)- 

(c) In order to deduce explicit analytical expressions for the modes, we expand the plane 
waves over vector spherical functions by using Eqs. ( |B6| ) . We then integrate the 
products of Wigner D-functions over the angles (pk and 6k by using the orthogonality 
condition ( [A1CI| ). 

(d) Finally, the modes are derived from the expressions for M^'^ by changing the 
Bessel functions, jj{p), to the Hankel functions, h^p{p). 

Note that, if a linear combination of Bessel functions, jj{p), represents a solution of linear 
homogeneous differential equations (Maxwell equations in our case), then the corresponding 
linear combination of Hankel functions generates another solution. This remark justiffes the 
last step in the procedure described above. 

The crucial point is that this inverse procedure can be generalised to a uniformly 
anisotropic medium. Thus it can be applied to superpositions of plane waves representing 
solutions of the Maxwell's equations in the uniformly anisotropic layer, yielding expressions 
for the modes used in the generalised T- matrix ansatz (p6|). In the next subsection we 
perform this generalisation. 



2. Wave functions in an anisotropic medium 



We start with the expansion (0), and use it to derive formulae for generalised spherical 
harmonics in the anisotropic medium. The starting point is the well known result for plane 
waves 1211, |2^, |23|, El: 



E = (exp(ipek ■ f)^,;(k) [e^(k) + -^-^ sin 9k z] + exp(2pk ■ f)^j,(k) ey(k))^^ , (42a) 

A A A AAA 

H = - (exp(ipek ■ f) n-^ E^{k)ey(k) - exp{ipk- r)Ey(k)e^(k))^^ , (42b) 
A* 

1 + u 

where = nl{9k) = — 5-- and pe = ne{9k) p- 

1 + u cos^ Ok 

We now apply the procedure described at the end of the last section to the plane wave 



packets (^21). This gives a representation of the electromagnetic field in the form of the 
generalised T-matrix ansatz (p6D. Now, however, the modes no longer take the form (p7[), 
but rather the modified form: 
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P^iP^ r) = ^-^•(47r)-V^(2j + 1) ( Dltik) exp(zp,k ■ f ) [e.(k) + -— sin 9, z] - 

-^Dfj*(k)exp(zpk-f)e,(k))£, (43a) 
P51(P>f) = ^'^i^n)-WA2J + l){^Df2*i^) exp(zp,k ■ f) [e.(k) + ^ sin^z] + 

+ ^?2*(k)exp(zpk-f)e,(k))fc, (43b) 
q5.::)(p, f) = r^(47r)-Vvr(2j + 1) ( -^/^fi*(k) exp(^pek ■ f) e,(k) + 

+ Dj^i*(k)exp(2pk-f)e,.(k))£, (43c) 
Q5-1(p, r) = r^-(47r)-V^(2j + l)( /^^^ *(k) exp(.pek ■ f) e,(k) + 

+ tD\t{^) ewi^P k ■ f) e.(k) )^ , (43d) 

Typically, solving a scattering problem for a spherical particle requires modes to be 
expressed in terms of vector spherical harmonics. It thus turns out to be useful to write the 
wave functions (^) as expansions over vector spherical harmonics: 



fB j'>\m\ 13 j'>\m\ 

13 j'> \m\ 13 j'> \m\ 

where a G {m, e}, /5 G {m, e, o} and 

p^^::lip) = (Y/i*(f) ■ pj:^(p, ?)), , pj^'^^V) = (Y/i*(f) ■ pj^^p, , 

Q^^:^ip) = (Y,S*(f) • Qi:^(P> , ^^;™^(P) = (Y,^i*(f) ■ Q,t^(P, f)). • (45) 

Explicit formulae for the coefficient functions entering Eqs.. (|44b|) are given in Appendix 
together with some related comments. Evaluating these coefficients involves computing 
some products of Bessel spherical function and Wigner D-functions and integrating these 
expressions over 6k- Their numerical evaluation is relatively easy. 

The coefficients with j ^ j' describe angular momentum mixing. It can be shown that 
these terms go to zero in the absence of anisotropy, when u = and = 1. 

The modes introduced in this subsection can be used to derive expansions for incident 
plane wave in the anisotropic medium. In an isotropic material the corresponding expansions 
are given by Eqs. (^ and (|^). In a uniformly anisotropic medium these expansions take an 
exactly analogous form: 

p(mc) ^ (mc) pM/ t.) ^ P n^'-"-'''^ P^*"^ (n f) 
jm ^jrn jm Kri '- ) ^ ^jm jm\l i ) i 

jm = "im QjJ (p. r) + - a}„ Q;„;(p, r) . (46) 

P 

These equations can be compared to Eqs. (^. In the anisotropic case the spherical harmon- 
ics d^) are replaced by by the "quasi-spherical" modes (^). 

In order to prove this result we insert the coefficients from Eqs. (|^) into Eqs. (p|)- (|4ff5D . 



Then, after performing some rather straightforward algebraic manipulations, we obtain 
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E(-^)(k.„„ k) = ^ ^(2j- + l)/(47r) z/E^^)D^,(k_)Z}^^(k) , (47a) 

u=±l jm 

Ef'=)(k,.„ k) = -t 2-1/2 ^ ^(2j + l)/(47r) 4^"^)l^4.(k.nc)I^^*.(k) . (47b) 

v=±l jm 

We finally observe that the set {Dl^j^ij^ }jm form an orthogonal (see Eq. ( [A10| ) )and 
complete set of functions in the space of angular dependent functions. The sums on the 
right hand sides of Eqs. ( |47a| ) and ( [47b| ) can thus be written as angular 5-functions: 

^(2j + l)/(47r) D^,(k,„,)D^*,(k) = 5(k - k,„,) . (48) 

jm 

B. T-matrix: uniform anisotropy 

In Sec. [IV A| we solved for the T- matrix of a radially anisotropic layer. Computing the 



elements of T-matrix required the solution of a set of equations resulting from the boundary 
conditions (|29|) . The only difference in other cases is that that appropriate modes for the 
electromagnetic field inside the anisotropic layer must be used. For uniform anisotropy these 
modes are given by Eqs. (|4^) . Mathematically, Eqs. (|30|) and (|3lD , which describe the fields 
in radially anisotropic layer, are replaced by the following relations: 



where 



TjmKn 
tm) I ■ 
IjmV. 

\vfmir)) 



j'>\m\ 



(r) 



/ 0!j'm\ 



f^j'm 
yPj'mJ 



(49) 



R 



33 vmi 



( 



V 



(m,m) 
jj';m 



(P 



nijJi q 



1 ^(e,e) 



1 ^(m, 



^iPl ^ <ljj'-m 
\ (e,m) / \ 

\ Pjj';m{Pl) 



(Pl 
(Pl 



~ (m,m) 



(Pl) 

1 Hjji-miP'^/ 
(Pl) 



-1 ^(e,e) 

~ (e,m) 
Pjj'; m 



-1 



(Pi 



— 1 (m,e) 

(e,m) / \ 
<ljj';m{Pl) 

4™C^(Pl) 

— 1 (e,e) / ^ 



— 1 ~(m,e)/ \\ 
'^1^1 :Pii';rn(Pl)\ 

1jj'-m\Pl) 
~{m,m), N 
?jj';m iPlJ 

-1 ~(e,e) 



(Pi)/ 



(50) 



These changes affect the left hand sides of the system (|3^ ) which is modified in the following 
way 
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j' > \m\ 



Rj3';m 
2 



j' > |m 



Ctj 

V4- 



( Olj'm' 



(c) 
jm 



(sea) 
jm 



, (inc) 

+ a • 



/ [ji(P2)]2 \ 



V / 

/ [hf\p)]i \ 
n/fi[hf\p)][ 


V 

/ bi(p)]i \ 

n/^J'[jjip)]'l 





, ~(c) 






_bi(P2)]2 
V-^2 V2bj(P2)]2y 

/ 



(51a) 



(sea) 



J 



~ (ine) 

a ■ 

jm 



[hf\p)]i 



+ 



( 



\-Pln[h';\p)y 

\ 



bi(p)]i 

'^^/n[3j{p)]'i) 



(51b) 



By analogy with the procedure adopted in Sec. [IV A|, this system of equations can now 



be solved, in principle, to yield solutions for the quantities Tj"'.^ which occur in Eq. ([TT]). 
However, the crucial difficulty in this case is that the algebraic structure of the equations 
is complicated by the presence of angular momentum mixing. Thus, by contrast with the 



radially anisotropic layer considered in Sec. [IV A| , we are now unable to derive expressions 
for the elements of the T-matrix in closed form. The solution of this system of equations 
now requires numerical analysis. In future publications in this series, we shall discuss this 
problem in greater detail, and present some explicit results. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have developed a T-matrix approach which can describe light scattering 
by spherical scatterers containing optically anisotropic material arranged in an annular layer. 
We have confined our discussion to the two cases, which we have called, using natural 
language, radially and uniformly anisotropic systems. Just as in the related case when the 
scatterer itself is anisotropic, but the scattering material is optically isotropic, the presence 
of optical anisotropy affects the algebraic structure which underlies the T-matrix theory. 

From a mathematical point of view, the radial and uniform anisotropics present inter- 
esting but contrasting features. Here we draw the reader's attention to some of the most 
striking of these. 

The simplest case is the radially anisotropic layer. Here the scattering material is locally 
optically anisotropic, but because of the way that the anisotropy is arranged, the scatterer 
itself is spherically symmetric and remains a globally optically isotropic object. The result 
is that there is no angular momentum mixing, although we do have to introduce the new 
normal mode structure ( p6D within the anisotropic layer. 

We shall discuss quantitative results in detail in a later paper in this series. However, it 
is our purpose here to show briefly that the anisotropy effects can be important, even for 
this relatively simple case. As an example, we compare scattering by a radially anisotropic 
layer and a radially isotropic layer of the same dimensions, chosen in some some sense to be 
the best guess to an equivalent isotropic scatterer. 

Speciflcally, we consider scattering by a scatterer with a radially anisotropic layer as 
discussed in section 2, with the e_L within the layer matching e in the core and outside the 
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FIG. 2: Scattered intensity (see Eq. (p4[)) versus the scattering angle at kR2 = 1.5, kd = 4.0 
and ui = 0.25 for (a) effective isotropic scatterer with m^.s ~ 1.05127 (Qeff = 0.14391). and 
(b) radially anisotropic layer, n = f {Qmis = 0.1439). Insert at the upper right corner enlarges the 
backscattering tail. 

scatterer, and anisotropy coefficient u inside the layer. The equivalent effective scatterer 
has the same dimensions, possesses an isotropic layer of refractive index ries, so that ees is 
the dielectric constant and = Ues/n is the optical contrast. The refractive index nes 
is chosen so as to match scattering efficiencies of the anisotropic layer and of the effective 
scatterer. 

We compare the angular dependence of the scattering and the depolarisation factors in 
the two cases. The aggregate scattering is the same, by definition. But disaggregated, we 
get different contributions, both when looking at angles and when looking at polarisation 
shifts. In Fig. ^ we compare the angular dependence of the intensities. In the case we 
have considered, the effective scatterer gives a much more forward scattering signature, and 
the relatively tiny backward scattering contribution has a very different angular structure. 
Likewise, we see from Fig. § that the depolarisation factor (pSf ) as a function of the scattering 
angle, 9sca, is also very sensitive to the presence of anisotropy. 

Given these relatively large effects for what might be thought of as minor anisotropy, 
there is every reason to suppose that the influence of a uniform anisotropic layer will be 
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T 




sea sea 



FIG. 3: Depolarisation factor (see Eq. (p5|)) versus the scattering angle for both the radially 
anisotropic layer and its effective isotropic scatterer at kR2 = 1.5, ui = 0.25 and (a) kd = 1.0 
(meff « 1.04512, Qeff = 0.010724, Qrad = 0.010729), (b) kd = 4.0 {m^s ~ 1-05127, Qcfr = 0.14391, 
Q™d = 0.1439). 



even more profound. In this case the light-scattering problem is not exactly soluble. The 
key point is that the exact solutions for uniformly anisotropic medium are known as plane 
waves, whereas the spherical shape of the particle requires using some kind of spherical 
modes. 

In the case of the uniform anisotropy we have found it necessary to examine the relation 
between spherical harmonic expansions and plane wave solutions of Maxwell's equations. 
We have found that, by choosing the appropriate basis in k-space, we can define 'quasi- 
spherical' normal modes. These modes are exact solutions of Maxwell's equations and as 
such mix different angular momentum. However, in the limit of zero anisotropy, these modes 
tend to familiar spherical modes. More importantly, these quasi-spherical modes turn out 
to be relatively easily accessible computationally. Thus, there is every reason to suppose 
that the strategy can be adopted in rather more complicated situations. 

One such problem is the light scattering problem for a Faraday-active sphere. This 
problem has been treated using perturbation theory in |^ to explain the origin of magneto- 
transverse light diffusion known as the "photonic Hall effect" p6| , p7| . 

We now try to place this problem in a more general physical context. We were first 
motivated by the technological problem of describing light transmission through media with 
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liquid crystalline inclusions, and the inverse problem in which the matrix is liquid crystalline 
but the scatterers are isotropic. There is considerable current interest in such materials for 
optical applications and displays. The complete problem of light transmission through such 
materials not only involves the single scattering processes discussed in this paper, but also 
more general multiple scattering processes. 

The T-matrix formalism is a natural language within which to discuss such problems, 
beginning with single scattering theories of the type discussed in this paper, one can in 
principle construct an effective medium theory using, for example, the coherent potential 
approximation (CPA) or coated CPA ||2^, These theories determine effective optical 
characteristics of the medium from the condition that the scattering cross section is minimal 
or equal to zero on average. Since this requires averaging over director orientations, it is 
important to use basis functions with well defined transformation properties under rotations. 

This paper is the first in series of papers designed to improve understanding of scattering 
by optically anisotropic bodies, both singly and as components within complex anisotropic 
media. In the next papers in this series, we shall carry out detailed calculations using the 
theory presented in this paper, and compare the results with results derived using the simpler 
Rayleigh-Gans (RG) and Van de Hulst or Anomalous Diffraction Approximations (ADA). 
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APPENDIX A: VECTOR SPHERICAL HARMONICS 

In this appendix we introduce notations and definitions used throughout the paper. In 
addition, we relate the vector spherical harmonics and Wigner Z)-functions. 
Let us define the vectors 



e±i(f) = t('.?±2<^)/v^, eo(f) = f, (Al) 

where (p = (— sin 0, cos0, 0), 'd = (cos cos 0, cos 6' sin 0, — sin 6') are the unit vectors tan- 
gential to the sphere; and 6 are Euler angles of the unit vector f. These vectors can be 
expressed in terms of the vectors of spherical basis, eo = z, e±i = =i=(x ± iy) / ^/2, (x, y and 
z are the unit vectors directed along the corresponding coordinate axes) as follows: 



M=-l 



u e {±1,0} 



(A2) 



= exp{—ifi(f)) d^j^{6) is the Wigner D-function [|T^, The 
following properties of the vectors (|A2D are easy to verify 



where D^,(f) ^ D^,(0, 



(e:(f) . e,(f)) = S, 



(A3) 
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The vector spherical functions Y^-^ from Eqs. (j^) are expressed in terms of the vector 



spherical harmonics Yijm [0 defined by 



Ijm 



'^ly^-uunJlm^u ® e^,, l=j + 6, 6 e {±1, 0} 



(A4) 



u=-l 



where Yim = ■\/(2/ + 1)/ (47r)D^*Q is the spherical function |]T9|] and C4t_Vjm denotes the 
Clebsch-Gordon (Wigner) coefficient, 



X — 0>j X j+lj m -r X m, x-^ — x^j 



■ jm 

y(°) 

jm 



(A5a) 

^ [j7(2j + c, ^ [(j + l)/(2j + . (A5b) 



Eqs. and the equality [[18 



k\k2k mk / j m\m2m^m\k\ 



32 

"12^2 ' 



(A6) 



mim2 



give 



2j + l 
47r 



1/2 



(A7) 



where Cj^' = [{21 + l)/(2j + l)Y/^Cl^li , so that the non-vanishing values of C^^ are 



V2C^' = -u, V2C^7'^' = = 9, v^Cit'' = C^~'' = sj . (A8) 

From Eqs. ( |A5| ) , ( |A7| ) and ( |A^ ) we express the vector spherical harmonics in terms of the 
Wigner Z?-functions as follows 



Yjr(f) = [(2j + l)/87r]V2 {Z}^:_,(f) e_i(f) - D^„:i(f) e+i(f)} 
Yfi(f) = [(2j + l)/47r]V2z)^:o(f)eo(f). 



(A9a) 
(A9b) 
(A9c) 



Note that the D-functions meet the following orthogonality relations [1^, ^ 



47r 



2j + l 



(AlO) 



where (/)? = / d0 / sin^d^/. The orthogonality condition (|A10|) and Eqs. (|A9a|) 
Jo Jo 

( |A9c| ) show that a set of vector spherical harmonics is orthonormal: 



(All) 
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APPENDIX B: RAYLEIGH EXPANSIONS FOR VECTOR PLANE WAVES 



In this Appendix we comment on the vector version of the well known Rayleigh expansion 
(see, for example, p[): 



oo I 



exp(ik ■ r) = 47r ^ ^ i'ji{p) Yi^{r) Yil{k), p = kr . 



(Bl) 



/=0 m=-l 



Let us consider plane wave with the wave vector k = fck and the polarisation vector E 
defined by its components, E^, in the basis e,^(k) (see Eq. 



e.(k)= 5^D;,(k)e^, e±i(k) = T(e.(k)±ze,(k))/V2, E=^E.e,(k), (B2) 



Dj{^{k) = D^(j)k,ek) = exp{-iij(f)k) d^Ok) is the Wigner D-function 1T|, ^ and 0^, 
6k are the azimuthal and polar angles of the unit vector k. The vectors ey(k) = c^^ = 
(— sin0^., cos0fc, 0) and ea,(k) = 1?^ = (cos6'a; cos^^, cos6'fc sin 0^, — sin6'fc) are perpendicular 
to k. 

From Eq. (|B1|), definition of the vector spherical functions ( |X5| ) and the equality ( |A6|) it 
is not difficult to derive the following relation 



.(k)exp(zk-r) = ^[27r(2j + l)]i/^D^,(k) J^^' j/(p) Y,,^(f) 



(B3) 



where Cj-^ is defined in Eq. (|A8|) . The sum in square brackets on the right hand side of 
Eq. ( [B3D can be simplified by making use of Eq. ( |X8| ) and the recursion relations [0: 



dj- 

jj+i + jj-i = (2j + 1)P"^ jj, Jj-i - jj+i = 2-7^ + jj . 



dp 



The final result for transverse waves with z/ = ±1 is 



(B4) 



(k) exp(z k ■ r) = J^^^-' [2n{2j + 1)]'/' DU^) ■ 
It can be written in the following form 



e.(k)exp(^k-r) = 5^z^[7r(2j + l)]V2 

j,m 

e,(k)exp(^k-r) = 5^zM^(2j + l)]'/' 



^fi(k)M;-^(p,f)-zD|.^j(k)M(:^(p,f) , (B6a) 



jm \ / jm \ 



(B5) 



iDj:i(k)M(S(p,f) + 4£(k)M(:)(p,f) , (B6b) 
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where the mo des m5™\ Mfl are defined by Eq. (|5a| ) and the functions Dj^J^, Dj^ are 
expressed in terms of Wigner D-functions in Eqs. ( P8a| )-( ^8BD . 



Let us consider the case when the polarisation vector is directed along the z axis. We 
derive the formulae for the following matrix elements: ( Y^-^^* ■ zexp(zk ■ r) In order to 



do it, notice that the definition ( [A^ ) immediately gives the following relation: 



Ijm ' 



(B7) 



Then we substitute Eq. ( pTj) into the expansion multiplied by z and express the vector 
functions Y^j^ in terms of the vector spherical harmonics Yj^\ The result reads 



r(m) 
■ jm 

^(e)* 
- jm 



Y\-J* . z exp(zk • r) ), = ^^4n{2J + l)r/'CillDU^) j,{p) , 

■ z exp(z k ■ r) ), = z^+^[47r] V2 [ (2j + 3)^/^ C^+l'^' D^^o\^) jMp) ' 
- (2j - iy/'c,Ci~,'J,^D^,\k)j,^,{p) 
( . z exp(z k ■ r) ), = z^+'[Anf' \ (2j + 3)^/^ c, C^+U^' Di^,\k) j,^,{p) + 



+ (2j - 1) V2 C^oU^' ^^o'(k) J,-i(p) 



fB8a) 



(B8b) 



(BSc) 



with the Wigner coefficients given by |[T^, [T9| 

1/2 



m m 



m 



,j(2j 



m 



VJU + 1) 



L(j + l)(2j + 3), 



1/2 



(B9) 



In conclusion, note that in case, when a series on the right hand side of Eq. ( P5| ) or 
Eqs. (P8| ) represent a solution of the Maxwell equations (^, we can obtain another solution 
by replacing jj{p) with hf\p). 



APPENDIX C: COEFFICIENT FUNCTIONS 



By definition, the coefficient functions that enter the expansions p^ ) are the matrix 
elements (|^). In order to deduce the corresponding formulae we need to substitute the 
expansions for plane waves Eqs. (BB) into Eqs. and make use of the matrix elements 
for the plane wave polarised along the z axis given by Eqs. ( |B8a| )- ( [B8c| ). Note that Eqs. 
combined with the orthogonality condition for the vector spherical harmonics ( |A11| ) imme- 
diately give the following relations: 



(Yj:)*(f)-M(fi,(p,f))j 

'Y^{i).M%{p,r)),^ 



■ jm 



^m,a ^jj' ^mm' jj (p) 



(CI) 

(C2) 



Below we write the resulting expressions for the matrix elements that correspond to 
transverse components of the wave functions (^31) . 
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p^7-::^ip) = ^ ( j.(Pe) + c^j^^:^ j,(p) + a4:;'^(p) , (C4) 

pI;;^ (p) = ( ^i^^l ^j.(pe) + rfj^^^i ^j.(p) + ^P^r;:^iP) , (cs) 

p1;';1(p) = A^.' ( mpe) + d^;^i dj,{p) + ^p^:i{p) , (ce) 

^'r^(p) = ( ^j;::;! j.(pe) + d^^^tj^p) , (C7) 

(p) = ^ A^.' ( J. (pe) + rfj;;^ J. (P) , (C8) 

(P) = ( df^S^ Dj.ip^) n-' + Dj,{p) )e, , (C9) 

^;1(P) = ^N,,{d^:^Dj,{p:) n-' + rfj^;:;; Dj,(p) , (CIO) 



where ^ '-^[{2j + l)(2j' + 1)]V2 and dj;;^) ^ rfi:Vfc)4?l(^'^)' ^ 
The terms Apj^'.(2(p) given by 



Vj(j + l) 



A P^-tl (P) = ^^';1(P) , (C12) 



2jTl 



ApSC^ (P) = ^ b-((j + 1)^ - m^)l{3 + ir'r^'^,%Jp) 



- [{j + i){f - "^')/jr/v]!f),^(p) \ , (C13) 

2j + J 

- [(j + - m')/j]'/'r^%Jp) \ , (CM) 



^P^P:LiP) = ^ + 1)^ - + l)]^/^r(:'5,^(p) 



where 



Expressions for PjJ>'.^{p) and qjp.mip) can be derived from Eqs. (|C3|) -( |CT5D by replacing 
spherical Bessel functions of the first kind jj{x) with spherical Bessel functions of the third 

kind h^^\x). 

The coefficient functions Pjj'/'-m ^^"^ ^jf-m diagonal, oc da^p^jj', in the limit of weak 
anisotropy, n ^ 0. It can be seen from orthogonality condition ([A10|) that provides the 
orthogonality conditions for d-functions defined by Eqs. ( P8a|) - (|38b| ) in the following form 

\ dji-.rn + djij-rn/Sk = ' 



{df-'L + d^;t)e, = ^. (C16) 
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